Extreme points of the set of density matrices with positive partial transpose 
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We present a necessary and sufficient condition for a finite dimensional density matrix to be an extreme point 
of the convex set of density matrices with positive partial transpose with respect to a subsystem. We also give 
an algorithm for finding such extreme points and illustrate this by some examples. 



The density matrices describing states of a bipartite quan- 
tum system are uniquely classified as being either separable 
or entangled. However, it may be a non-trivial task in prac- 
tice to classify a given density matrix. One simple method 
which can sometimes prove a density matrix to be entangled, 
and which is especially useful in systems with Hilbert spaces 
of low dimensions, is due to A. Peres fl]]. Peres noted that a 
separable density matrix must remain positive when a partial 
transposition is performed with respect to one of the subsys- 
tems. In general the set of density matrices that have positive 
partial transpose, called PPT states for short, is larger than 
the set of separable matrices, but in low dimensions the dif- 
ference between these two sets is small. In particular, for a 
system composed of one subsystem of dimension 2 and one 
subsystem of dimension 2 or 3, the two sets are known to be 
identical HH. 

A particular approach to the question of separability is to 
characterize geometrically the different sets of matrices as 
subsets of the real Hilbert space of hermitian matrices 0]. 
Thus, the set of all density matrices is a compact convex set 
T>, with the set of separable density matrices, S, and the set of 
PPT matrices, V, which we will call the Peres set, as convex 
subsets. The Peres criterion states that S C V C V. We have 
presented elsewhere a numerical method for deciding whether 
a density matrix is included in S by computing the distance to 
the closest separable density matrix J^Hg]. 

There are two complementary ways to characterize a con- 
vex set. One way is to specify its extreme points. Thus, every 
point in a compact convex set of finite dimension d has an 
expansion as a convex combination of d + 1 extreme points, 
while an extreme point cannot be written as a convex com- 
bination of other points in the set. The other way is to iden- 
tify the set by conditions, typically algebraic equations and 
inequalities, that define whether or not a given point belongs 
to the set. For example, the convex set V may be defined by 
its extreme points, which are the pure states (one dimensional 
projections), or equivalently by the conditions of hermiticity, 
positivity and unit trace. In general there is no simple relation 
between these two descriptions of a convex set. 

It is interesting to note that the convex sets S and V have 



simple charcterizations in the complementary ways just de- 
scribed. Thus, the extreme points of S are well known, they 
are the pure product states of the two subsystems, but there is 
no (known) effective test for membership of S. In contrast, 
testing for membership of V is easy, because V = V n T> p , 
where the superscript P denotes partial transposition, but the 
extreme points of V are not completely known. 

The purpose of the present paper is to study the difference 
between the two sets S and V by studying the extreme points 
of V. The extreme points of S are also extreme points of V, 
because they are extreme points of V and S C V C T>. But in 
general V is larger than S and has additional extreme points, 
which are unknown and which make the whole difference be- 
tween the two sets. 

The extreme points of V that are not pure product states are 
interesting also as examples of entangled PPT states. In fact, 
because S C V, a separable extreme point of V must be an 
extreme point of S and hence a pure product state. It is easy 
to verify that any pure state which is not a product state, is not 
in V and is therefore entangled. Thus, every extreme point 
of V which is not a pure product state, is entangled and has 
rank higher than one. A stronger lower bound on the ranks of 
entangled PPT states is actually known J3]. 

We will in the following specify a criterion for uniquely 
identifying the extreme points of V and we will describe an al- 
gorithm for finding such points. The method is demonstrated 
by some examples. 

To be more specific we consider a composite quantum sys- 
tem with Hilbert space TL = Ha ® H-B of finite dimension 
N = Na Nb, where A and B denote the two subsystems. 
Then V is the set of density matrices on H, 

p<^V^p = p\ p>0, Trp=l (1) 

A subset is the Peres set V, 

p G V <=> p £ V and p p e V (2) 

where p p is the partial transpose of p with respect to any one 
of the subsystems, say system B. A subset of V is S, the set 
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of separable density matrices, 

P eS^p = Y,Pkp£®Pk, Pk>0, J2 Pk = 1 (3) 

k k 

Thus, <S is the convex hull of the product states of the two 
subsystems, and its extreme points are the pure product states. 

To develop the method we need for finding the extreme 
points of V we will first apply it to T>, the full set of density 
matrices, with the pure states as extreme points. We recall 
some definitions and elementary facts. 

Every p G D is hermitian and has a spectral decomposition 

p = Y^Xi\il>i)(il>i\ (4) 

i 

with real eigenvalues A; and orthonormal eigenvectors \tpi). 
The positivity condition p > means that all Xi > 0, or 
equivalently that (ip\p\ijj) > for all with (ip\p\i/j) = if 
and only if p\tp) = 0. The orthogonal projection 

i,Ai>0 

projects onto the image (or range) of p, whereas 1 — P projects 
onto the kernel of p, denoted by ker p. 

If p is not an extreme point of I? it is a convex combination 

p = xp' + (1 - x)p" , < x < 1 (6) 

with p 1 , p" e I? and p' ^ p". The identity 

(V>|/#> = a: <VIP'I^> +(!-») (7) 

shows that p > when p' > and p" > 0, thus proving the 
convexity of T>. More interestingly, it shows that 

ker p = ker p' n ker p" (8) 

With P defined as in (0 we have therefore PpP = p, Pp'P = 
p', and Pp"P = p". 

When © holds, the matrix a = p' — p is hermitian and 
nonzero, and Tr a = 0, hence er has both positive and negative 
eigenvalues. Moreover, Pop = a. Define 

t(x) = p + xa (9) 

for x real. Since a has both positive and negative eigenval- 
ues, so has t(x) for large enough If p\ip) = then 
P\tp) = and a\f) = PoP\4>) = 0, hence T(x)\tp) = 
for all x. Therefore only the strictly positive eigenvalues of p 
can change when xa is added to p, and since they change con- 
tinuously with x, they remain positive for x in a finite interval 
about x = 0. 

We conclude that there exists an x\ < and an x^ > such 
that t(x) > for x\ < x < X2, and r(x) ^ for x < x\ 
or x > X2- At x = x\ or x = X2, t(x) has at least one zero 
eigenvalue more than p. 

We are now prepared to search systematically for extreme 
points of T>. Starting with an arbitrary pi 6 T> we define the 



projection Pi in the same way as we defined P from p in (0. 
If we can find a hermitian matrix a solving the equation 

PicrPi = a (10) 

we define 

(7i=o-- (Trer)pi (11) 

in order to have P\0\P\ = <j\ and Tr o\ = 0. Clearly a = p\ 
is a solution of ( TTOb . If this is the only solution, then only a\ = 
is possible, and it follows from the above discussion that 
pi is an extreme point. The number of linearly independent 
solutions of ( [Tol l is n^, where n\ = Tr Pi is the rank of pi. 
Hence, pi is an extreme point if and only if ri\ = 1 so that it 
is a pure state. 

If pi is not an extreme point, then we can find o\ ^ and 
define 

n(x)=pi+xai (12) 

We increase (or decrease) x from x = until it first hap- 
pens that n (x) gets one or more additional zero eigenvalues 
as compared to pi. We will know if we go too far, because 
then n (x) will get negative eigenvalues. We choose p2 as the 
limiting t\ (x) determined in this way. By construction, pj has 
lower rank than p\ . 

We repeat the whole procedure with p 2 in place of p x , and if 
P2 is not extreme we will find a p3 of lower rank. Continuing 
in the same way, we must end up at an extreme point pk, with 
K < N, since we get a decreasing sequence of projections, 

I 2 Pi DP> D ••■ dPk (13) 

of decreasing ranks N > m > ri2 > . . . > %• = 1. 

We may understand the above construction geometrically. 
In fact, each projection Pfc defines a convex subset P^PPfe 
of T>, with pk as an interior point. This subset consists of all 
density matrices on the rifc dimensional Hilbert space PkH, 
and if < N it is a flat face of the boundary of T>. 

It is straightforward to adapt the above method and use it 
to search for extreme points of P. Thus, we consider an ini- 
tial matrix pi G P, characterized by two integers (ni,mi), 
the ranks of pi and of the partial transpose pf , since pi 6 P 
means that pi G T> and pf G V. We denote by Pi the pro- 
jection on the image of p%, as before, and we introduce Qi as 
the projection on the image of pf . We have to solve the two 
equations 

P 1 aP 1 =a, Qi<j p Q 1 =<j p (14) 

To understand these equations it may help to think of a as a 
vector in the iV 2 dimensional real Hilbert space A4 of N x N 
hermitian matrices with the scalar product 

(A, B) = (B, A) = Tr(AB) = ^ A* B tJ (15) 

i,3 
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If L is a linear transformation on Ai, its transpose L T is de- 
fined by the identity (A, L T B) = (LA, B), and L is symmet- 
ric if L T = L. Partial transposition of A G Ai permutes the 
matrix elements of A and is a linear transformation TIA = 
A p . It is its own inverse, and is an orthogonal transformation 
(it preserves the scalar product), hence II = II -1 = II T . It 
also preserves the trace, Tr(IIA) = Tr A. The projections Pi 
and Qi on H define projections Pi and Qi on Ai by 

V X A = P X AP X , QiA = QxAQx (16) 

These are both orthogonal projections: Pj 2 = Pi, = Pi, 
Qi 2 = Qi, and Qf = Qi 

In this language the equations ( fT4b may be written as 

Pl(T = cr, Qxa = cr (17) 

with Qi = EEQiII. Note that Qi is also an orthogonal pro- 
jection: Qi 2 = Qi and QJ = Qi. These two equations are 
equivalent to the single equation 

PlQlPl < T = ( T (18) 

or equivalently QiPiQicr = a. They restrict the hermitian 
matrix a to the intersection between the two subspaces of Ai 
defined by the projections Pi and Qi. We shall denote by 
Bi the projection on this subspace, which is spanned by the 
eigenvectors with eigenvalue 1 of P1Q1P1. Because the lat- 
ter is a symmetric linear transformation on Ai it has a com- 
plete set of orthonormal real eigenvectors and eigenvalues. All 
its eigenvalues lie between and 1 . We may diagonalize it in 
order to find its eigenvectors with eigenvalue 1. 

Having found a as a solution of (fT8l we define a\ as in (fTTb . 
If a = pi and ax = is the only solution, then p\ is an 
extreme point of V . If we can find a\ ^ 0, then we define 
Tx (x) as in (1121 . and increase (or decrease) x from x = until 
we reach the first value of x where either Tx(x) or (ti(x)) p 
has at least one new zero eigenvalue. This special t\ (x) we 
take as p2- By construction, when n-2 is the rank of pi and 
77i2 the rank of p P , we have 712 < nx, Ti2 < mi, and either 

772 < TlX Or 7772 < fTll- 

Next, we check whether P2 is an extreme point of V, in the 
same way as with pi. If P2 is not extreme, then we can find a 
third candidat p$, and so on. We will reach an extreme point 
Pk in a finite number of iterations, because the sum of ranks, 
n k + m k> decreases in each iteration. 

The geometrical interpretation of this iteration scheme is 
that the density matrices p <E V which satisfy the condition 
Bfe p = p define a convex subset of V having pk as an interior 
point. This subset is either the whole of V, or a flat face of 
the boundary of V, which is the intersection of a flat face of 
V and a flat face of T> p . 

The construction discussed above defines an algorithm for 
finding extreme points of V, and gives at the same time a nec- 
essary and sufficient condition for a density matrix in V to be 
an extreme point. Let us restate this condition: 

A density matrix p G V is an extreme point of V if and 
only if the projection P which projects on the image of p and 



the projection Q which projects on the image of p p , define a 
combined projection B of rank 1 in the real Hilbert space Ai 
of hermitian matrices. 

The algorithm deserves some further comments. In itera- 
tion k the projections Pk, Qk and B^ are uniquely defined by 
the density matrix pk, but the matrix is usually not unique. 
In the applications discussed below we have simply chosen 
CTfc randomly. Clearly, more systematic choices are possible if 
one wants to search for special types of extreme points. 

Another comment concerns the ranks (77, 777) of a density 
matrix p and its partial transpose p p if it is an extreme point. 
We can derive an upper limit on these ranks. The rank of the 
projection P is tt 2 and the rank of Q is ttt, 2 , thus, the equations 
Per = <7 and Q<7 = a for the A 2 dimensional vector a repre- 
sent A 2 — 77 2 and A 2 — m 2 constraints, respectively. The total 
number of independent constraints is n c < 2A 2 — n 2 — m 2 , 
and the rank of B is A 2 — T7 C > 77 2 + m 2 — A 2 . For an extreme 
point the rank of B is 1, implying the inequality 

77 2 +m 2 <A 2 + l (19) 

We have used our algorithm in numerical studies. In one ap- 
proach we use as initial density matrix the maximally mixed 
state pi = 1/A and choose in iteration k a random direc- 
tion in the subspace B^AL We list in Table 1 all the ranks 
of extreme points found in this way in various dimensions 
A = NaNb- We do not distinguish between ranks (n,m) 
and (to, 77) since there is full symmetry between p and p p . 
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TABLE I: Typical ranks (n, m) for extreme points 

We find only solutions of maximal rank, in the sense that 
increasing either norm will violate the inequality (fT~9b . Max- 
imal rank means that the constraints given by the equations 
Per = (7 and Qcr = a are mostly independent. Further- 
more, we find only the most symmetric ranks, in the sense 
that m w 77,. For example, in the 4x4 system we find ranks 
(11, 11) and (10, 12), but not (9, 13), (7, 14) or (5, 15), which 
are also maximal. 

The 3x3 system we have examined further in the following 
way. For one specific sequence pi , P2 , . . . , pi< with px ex- 
treme, we repeat the final step, keeping pk-i fixed but choos- 
ing different directions in the subspace Ba-_iA4. We find 
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FIG. 1: Section through the boundary of V, in 3 x 3 dimensions. A 
closed curve of extreme points surrounds a region of entangled PPT 
matrices of rank (7, 6). The curve has two parts, characterized by 
ranks (7, 5) and (6, 6), joining at two rank (6, 5) extreme points. 

that every direction points directly towards an extreme point. 
Thus, pk-1 is an interior point of a flat face of the bound- 
ary of V bounded by a hypersurface of extreme points. Fig. 1 
shows a two dimensional section through this flat face. 

This shows that extreme points of non-maximal rank do ex- 
ist, and the fact that we do not find them in random searches 
just indicates that they define subsets of lower dimension than 
the extreme points of maximal rank (n, m) with m k, n. Note 
that an extreme point which is not a pure product state cannot 
have arbitrarily low rank, since in |7|] there is a proof that all 
PPT matrices of rank less than A*o = min{iVA, Nb} are sep- 
arable. This implies a lower limit of (No, No) for the ranks 
of an extreme point of V which is not a pure product state. It 
is not known whether there exist entangled PPT states of the 
minimal rank Nq. 

Several of the low rank entangled PPT states that are known 
turn out, in our test, to be extreme points. Examples in 3 x 
3 dimensions include the unextendible product basis state of 
rank (4,4) H; another state of rank (4,4) U; and explicit 
examples of rank (5, 5) and (6, 6) states II 1 011 . 

The entangled PPT states first discovered 13], in 3 x 3 di- 
mensions to be specific, are not extreme points of V, but on the 
flat face defined by the corresponding projection B they seem 
to be completely surrounded by extreme points. Figure 2 is a 
two dimensional section chosen so as to show one such state 
(with parameter value a = 0.42, called here the "Horodecki 
state"), as a convex combination of two extreme points of V . 
We would expect a two dimensional section through two ex- 
treme points of V to show maximum difference between the 
sets S and V . Thus, this plot illustrates the fact that the differ- 
ence is indeed very small in 3 x 3 dimensions. 

In conclusion, the method discussed has the potential of 
producing a clearer picture of the difference between the two 
sets S and V and thereby the set of states with bound entan- 
glement. We intend to follow up the work presented in this 



paper by other numerical studies of composite systems of low 
Hilbert space dimensions. 




FIG. 2: Section through the set of density matrices in 3 x 3 dimen- 
sions. The star in the middle of a straight line is the "Horodecki 
state" (see text). It is a convex combination of two extreme points, 
one of which is plotted as a star. The maximally mixed state 1/N 
is the star close to the center. The separable matrices lie in the large 
(red) region marked S, the entangled PPT states in the (purple) re- 
gion marked B, and the entangled non-PPT states in the two (blue) 
regions marked E. The lines are solutions of the equations det p = 
(blue) and det p p — (red lines). 
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